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ABSTRACT 

We present results from spectral modeling of three black hole X-ray binaries: LMC X-3, GRO 
J1655-40, and XTE J1550-564. Using a sample of disk dominated observations, we fit the data with a 
range of spectral models that includes a simple multitemperature blackbody (DISKBB), a relativistic 
accretion disk model based on color-corrected blackbodies (KERRBB), and a relativistic model based 
on non-LTE atmosphere models within an a prescription (BHSPEC). BHSPEC provides the best fit 
for a BeppoSAX observation of LMC X-3, which has the broadest energy coverage of our sample. 
It also provides the best fit for multiple epochs of Rossi X-ray Timing Explorer (RXTE) data in 
this source, except at the very highest luminosity {L/L-Edd ^ 0.7), where additional physics must be 
coming into play. BHSPEC is also the best-fit model for multi-epoch RXTE observations of GRO 
J1655-40 and XTE J1550-564, although the best-fit inclination of the inner disk differs from the binary 
inclination. All our fits prefer a = 0.01 to a = 0.1, in apparent disagreement with the large stresses 
inferred from the rapid rise times observed in outbursts of these two sources. In all three sources 
our fits imply moderate black hole spins (a* ~ 0.1 — 0.8), but this is sensitive to the reliability of 
independent measurements of these system parameters and to the physical assumptions which underly 
our spectral models. 

Subject headings: accretion, accretion disks — black hole physics — X-rays:binaries — Stars: Individ- 
ual(GRO J1655-40, LMC X-3, XTE J1550-564) 



1. INTRODUCTION 

Black hole X-ray binaries (BHBs) are known to occupy 
distinct spectral states which can be characterized by the 
relative contribution of thermal and non-thermal emis- 
sion components (e.g. McClintock & Remillard 2004). 
The most well-understood of these states is the thermal 
dominant (or high/soft) state. Here most of the flux is 
in the thermal component, which is generally assumed 
to be emission from a radiatively efficient, geometrically 
thin accretion disk (Shakura & Sunyaev 1973). The disk 
is believed to extend deep within the gravitational field 
of the black hole, and this makes spectral modeling of 
this state an important probe of both the physics of rel- 
ativistic accretion disks and the properties of black holes. 

In standard treatments of black hole accretion disks, 
the emitting matter extends down to the innermost sta- 
ble circular orbit (ISCO), which is determined by the 
mass and spin of the black hole. Such treatments usu- 
ally assume a "no-torque" inner boundary condition at 
this radius (e.g. Novikov & Thorne 1973), but magnetic 
fields may in fact exert such torques (Gammie 1999; Kro- 
lik 1999; Hawley & Krolik 2002), increasing the radiative 
efficiency of the disk. The structure and emission of these 
disks are therefore sensitive to the mass and spin of the 
black hole as well as any torque which may be present. 

This sensitivity makes accretion disk spectral modeling 
a potential way to measure or constrain black hole spin 
(e.g. Ebisawa et al. 1993; Zhang et al. 1997; Shafee et 
al. 2006; Middleton et al. 2006). This method requires 
a model of the radial profile of effective temperature in 
the gravitational field of the black hole, calculation of 

^ Department of Physics, University of California, Santa Bar- 
bara, CA 93106 

^ Department of Physics, University of Durham, South Road, 
Durham, DHl 3LE, UK 



the relativistic transfer function from the disk surface to 
an observer at infinity (Cunningham 1975), and spectral 
modeling of the surface emission in the local rest frame 
of the disk. Most implementations of this method ap- 
proximate one or all of these components. One common 
approximation is to assume that the disk surface emis- 
sion is a blackbody or, more generally, a color-corrected 
blackbody, 

= r^B,{fT,s), (1) 

where is the specific intensity. Toff is the effective tem- 
perature, Bi, is the Planck function, and / is the spectral 
hardening factor (or color-correction) , typically assumed 
to be around 1.7 (Shimura & Takahara 1995). One of 
the most sophisticated models of this type is the KER- 
RBB model (Li et al. 2005) for Xspec (Arnaud 1996). It 
accounts for the relativistic effects on the disk effective 
temperature profile and the relativistic transfer function. 
Potential difhculties with the color-corrected blackbody 
approximation exist. The local spectrum may not be well 
approximated by an isotropic, color-corrected blackbody 
due to limb darkening and frequency dependent absorp- 
tion opacities. Even if it can, one must still specify /. It 
has been suggested that / is a relatively strong function 
of accretion rate or of the fraction of energy emitted in 
a corona (Merloni et al. 2000). 

Relativistic models (Davis et al. 2005) now exist which 
calculate the non-LTE vertical disk structure and radia- 
tive transfer self-consistently using the TLUSTY stellar 
atmospheres code (Hubeny & Lanz 1995). The relativis- 
tic effects on photon geodesies are accounted for with 
ray tracing methods (Agol 1997). These spectral mod- 
els have now been implemented in an Xspec table model 
(BHSPEC; Davis & Hubeny, 2006). By calculating val- 
ues of / appropriate for use with KERRBB, these mod- 
els have already been used to estimate black hole spins 
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from KERRBB fits to two BHBs: GRO Jf 655-40 and 4U 
1543-47 (Shafee et al. 2006). 

In this work, we circumvent the color-corrected black- 
body approximation entirely by fitting the BHSPEC 
model directly to BHB observations. Since the model 
does not include irradiation of the disk surface, we fo- 
cus our efforts on thermal dominant state observations 
in which the non-thermal emission is a small fraction of 
the total flux. Fortunately, a sample of such observations 
made with the Rossi X-ray Timing Explorer (RXTE) al- 
ready exists (Gierlinski & Done 2004; hereafter GD04). 

BHSPEC still assumes that the disk emission extends 
only down to the ISCO and then effectively ceases due 
to rapid in-fall of matter interior to this radius. This as- 
sumption appears to be consistent with spectral model- 
ing of BHBs in the thermal dominant state in that the lu- 
minosity L is seen to scale roughly with the fourth power 
of the color temperature Tc in several different sources 
(Kubota et al. 2001; Gierlinski & Done 2004). This sug- 
gests that as sources vary by over an order of magnitude 
in luminosity, there is a roughly constant emitting area 
and thus a relatively constant inner radius to the disk. 
It is therefore very natural to associate such a stable in- 
ner radius with the ISCO of the black hole, though the 
'emission edge' need not coincide exactly with the ISCO 
(Krolik & Hawley 2002). 

Not all of these sources follow the L (x relation ex- 
actly, however. In several cases a relative hardening is 
seen with increasing L (see e.g. GD04; Kubota & Mak- 
ishima 2004; Shafee et al. 2006). A potential explana- 
tion for this hardening is that advection may be becom- 
ing increasingly important as these sources approach the 
Eddington limit (Kubota & Makishima 2004). Alterna- 
tively, this behavior is qualitatively consistent with the 
increased spectral hardening with accretion rate in the 
local disk atmospheres, ignoring advection (Davis et al. 
2005; Shafee et al. 2006). In the BHSPEC model, the 
precise nature of the hardening depends strongly on the 
variation in surface density with radius in the disk. Cur- 
rently, the surface density is determined by assuming the 
vertically averaged stress is proportional to the vertically 
averaged total pressure with a constant of proportional- 
ity a. However, more general stress prescriptions could 
be implemented in the future. Therefore, spectral mod- 
eling could potentially provide a constraint on the nature 
of the stresses in these systems. 

In this paper we fit these fuUy-relativistic, non-LTE 
accretion disk models to RXTE and BeppoSAX obser- 
vations of BHBs in the thermal dominant state. Our 
purpose is two-fold: we want to test the applicability 
of the spectral models to these observations, and hav- 
ing then found suitable representations of the data, we 
use them to infer black hole spins and infer properties of 
the stresses in these system. We review our spectral fit- 
ting and results in sectional discuss the implications of 
these results in section|3| and summarize our conclusions 
in section 0] In the Appendix we develop a simplified 
model to understand and motivate the variation of the 
BHSPEC spectra with accretion rate. 

2. SPECTRAL MODELING 

Several properties of BHBs make them particularly 
well suited both for testing accretion disk theory and for 
measuring the unknown black hole properties of interest. 
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Fig. 1. — The unfolded spectrum for the BeppoSAX observation 
of LMC X-3 using the best fit BHSPEC models for i = 67° , D = 
52 kpc, and a = 0.01. The total model component (green, solid 
curve), BHSPEC (red, long-dashed curve) and COMPTT (violet, 
short-dashed curve) are plotted. The unabsorbed BHSPEC model 
(orange, solid curve) is also shown. 



There are several sources for which precise, independent 
measurements of the mass of the primary and the binary 
inclination are available from light curve modeling of the 
secondary star (e.g. Orosz & Bailyn 1997). Typically, 
the distance to BHBs are known with less precision, but 
there are exceptions. As can be seen in Tabled the dis- 
tances to LMC X-3 and GRO J1655-40 are both claimed 
to be known to better than 10%. BHBs have an ad- 
vantage over most active galactic nuclei because of their 
relatively short timescales for large changes in the bolo- 
metric luminosity, which allows the same source to be 
observed in the same spectral state at appreciably dif- 
ferent accretion rates. The variation of the BHSPEC 
model with accretion rate provides a precise quantita- 
tive prediction which is sensitive to the assumed stress 
prescription. Therefore, simultaneous fitting of multiple 
observations of the same source at different epochs pro- 
vides a much more powerful constraint and potentially 
more information than a single fit to a single epoch. 

We compare three accretion disk models for the soft, 
thermal emission: the multicolor disk model (DISKBB in 
Xspec, Mitsuda et al. 1984), KERRBB, and BHSPEC. 
DISKBB is the most commonly used spectral model, but 
it neglects relativity Both KERRBB and BHSPEC in- 
clude relativistic effects. In addition, BHSPEC includes 
atmosphere physics which we are interested in testing in 
this paper. In particular, we are interested in examining 
whether BHSPEC with a fixed a can explain the vari- 
ation in the spectral hardening with accretion rate. As 
a control, we compare these results with KERRBB at a 
fixed /. 

We also model the neutral absorption along the line of 
sight, and though we focus on observations inferred to be 
disk dominated, we need an additional component to ac- 
count for the non-thermal emission which is present. It is 
widely believed that the coronal emission is due to inverse 
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Compton scattering of seed photons from the accretion 
disk. In this case, a power law will tend to overestimate 
the flux at low energies. Therefore, we prefer to ap- 
proximate the non-thermal emission with the COMPTT 
Comptonization model (Titarchuk 1994) in our spectral 
fits. We specify a disk geometry and fix the high energy 
cutoff in this model to 50 keV, a value high enough not 
to affect our fits. We also tie the seed photon tempera- 
ture to the DISKBB model temperature when DISKBB 
is used to model the soft emission. For KERRBB and 
BHSPEC we fix this parameter at the best fit DISKBB 
value. This leaves two free parameters - an optical depth 
and a normalization for each data set. 

DISKBB is a relatively simple model with only two pa- 
rameters: the temperature at the inner edge of the disk 
Tin and the model normahzation. BHSPEC and KER- 
RBB both share a number of model parameters which 
need to be specified or fit. The black hole spin a* = a/M 
and the accretion rate M are free parameters in all fits 
unless stated otherwise. For BHSPEC M is parameter- 
ized by £ = L/L-Edd where L^dd is the Eddington lu- 
minosity for completely ionized hydrogen. This value is 
converted to an accretion rate by assuming an efficiency 77 
which corresponds to the fraction of gravitational binding 
energy at infinity which is converted to radiation. Other 
parameters include black hole mass M, disk inclination 
i, and distance to the source D. For KERRBB D is an 
explicit parameter, but for BHSPEC D = 10/VN kpc 
where N is the model normalization. In some cases M, 
i, and D (or N) are fixed at the estimates given in Table 
nj Though M is always fixed, there are also cases where i 
and D are left as free parameters. The estimates for i in 
Tableware all estimates of the binary inclination. How- 
ever, there is no guarantee that the angular momentum 
of the black hole, and therefore the inner accretion disk, 
is aligned with the binary (Bardeen & Petterson 1975). If 
one assumes that the jet axis is aligned with the angular 
momentum vector of the black hole, then misalignments 
may be common (Maccarone 2002). XTE J1550-564 and 
GRO J1655-40 are among the sources for which misalign- 
ment can be inferred. Therefore, we also consider fits in 
which i and D are free parameters with i unconstrained 
and D allowed to vary within the confidence intervals in 
Table These models also allow for the presence of a 
torque on the inner disk which is parameterized by the 
increase in efficiency due to the torque relative to the 
efficiency of the untorqued disk A77/77. 

There are other parameters which are not shared be- 
tween BHSPEC and KERRBB. Two additional param- 
eters for BHSPEC are a and the metal abundance. In 
the BHSPEC model, the stress is given by 

TR^ = aP (2) 
where r^^ is the vertically averaged accretion stress and 
P is the vertically averaged total pressure. Usually, the 
metal abundances are fixed at the solar value, but we 
also consider fits with three times solar metallicity. The 
color correction parameter / must be chosen for KER- 
RBB, and we fix / = 1.7 unless stated otherwise. The / 
value which brings KERRBB into best agreement with 
BHSPEC is a function of ^ (Shafee et al. 2006), but this 
choice makes KERRBB roughly consistent over the range 
of £ we consider. We always fix the parameters rflag 
and Iflag so that the spectra are limb darkened and re- 
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Fig. 2.— The 66%, 90%, and 99% confidence contours in the 
a«-cosi plane for the best fit BHSPEC model (a = 0.01, i free) 
to the BeppoSAX LMC X-3 data. The vertical dashed lines mark 
uncertainty limits inferred for the binary inclination. 



processed emission from self-irradiation is ignored. Test 
cases suggest that these choices do not have a significant 
effect on the quality of fit or the inferred values for a*. 

2.1. Source Selection 

Our work is motivated in part by spectral fitting of 
BHBs performed by Gierlinski & Done (2004). They 
already provided a sample of sources with RXTE ob- 
servations in which a low fraction (under 15%) of the 
bolometric flux is inferred to be in the non-thermal com- 
ponent. We focus on three of these sources: LMC X-3, 
XTE J1550-564 (hereafter J1550), and GRO J1655-40 
(hereafter J1655). Each one ranges over nearly a decade 
or more in bolometric luminosity (see Figure 2 of Gier- 
linski & Done 2004). This makes them particularly well 
suited for constraining the spectral variation over a range 
of M . The properties of these sources are summarized 
in Tabled All have reasonably precise mass estimates, 
and the distances to LMC X-3 and J1655 both have rel- 
atively small uncertainties. The distance to J1550 is less 
well constrained, but we still include it in our sample 
because it spans the widest range of luminosities. 

One drawback of RXTE is its lack of soft X-ray cover- 
age. The thermal components of BHBs typically peak at 
photon energies near or below 1 keV, whereas the RXTE 
band extends down to only 3 keV. Since one of our 
goals is to test the applicability of the underlying ac- 
cretion disk model, we would like to cover as much of 
the SED as possible. Even if we use other observato- 
ries, we face the difficulty that most BHBs lie in or near 
the Galactic plane and are heavily absorbed by the inter- 
stellar medium along the line-of-sight. Therefore, we also 
examine a BeppoSAX observation of LMC X-3 which has 
a low line-of-sight absorption column. BeppoSAX is well- 
suited for this purpose as it covers a very broad range of 
photon energies extending down to a tenth of a keV, but 
lacks the effective area of XMM or Chandra and their 
corresponding pile-up problems for such bright sources. 

2.2. LMC X-3 
2.2.1. BeppoSAX Data 
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TABLE 1 

Source Descriptions 



Name 


Mass 


Distance 


Inclination 




Reference 




(Mo) 


(kpc) 


(deg) 


(10^2 cm-2) 




LMC X-3 


7(5-11) 


52(51.4 - 52.6) 


67(65 - 69) 


0.04 


1, 2, 3, 4, 5 


XTE J1550-564 


10(9.7 - 11.6) 


5.3(2.8 - 7.6) 


72(70.8 - 75.4) 


0.65 


6, 7 


GRO J1655-40 


7(6.8 - 7.2) 


3.2(3.0 - 3.4) 


70(64 - 71) 


0.8 


8, 9, 10, 11 



References. — (1) Soria et al. (2001); (2) Cowley et al. (1983); (3) di Benedetto (1997); (4) Kuiper et al. (1988); (5) Page et al. 
(2003); (6) Orosz et al. (2002); (7) Gierlinski & Done (2003); (8) Shahbaz et al. (1999); (9) Hjellming & Rupen (1995); (10) van der Hooft 
et al. (1998); (11) Gierlinski et al. (2001) 



The low energy coverage of BeppoSAX makes the fits 
particularly sensitive to our model for the neutral ab- 
sorption. We therefore include the line-of-sight absorp- 
tion column as a free parameter in our fits to the Bep- 
poSAX data. For our best fit BHSPEC model we find 
Nh = 5.92toiY X 10^0 cm-2. This is higher than values 
inferred from 21 cm absorption (3.2 x fO^'^ cm~^, Nowak 
et al. 200 f) or fits to the neutral oxygen edge in observa- 
tions with the Reflection Grating Spectrometer on-board 
XMM-Newton (3.8±0.8 x lO^o cm^^ Page et al. 2003), 
but lower than those found in previous modeling of the 
BeppoSAX data (7±1 x lO^o cm'^, Haardt et al. 200f). 

The unfolded spectrum of the BeppoSAX data is shown 
in FigureHl The best fit model for BHSPEC with i = 67° 
and a — O.Of which was used to generate the unfolded 
spectrum is also shown. The LMC X-3 spectrum is very 
disk dominated and the PDS provides little constraint 
due to the low count rate. Therefore, we fix the optical 
depth in the COMPTT model at r = 0.5, providing a 
flat (in vF^) spectrum typical of thermal dominant state 
observations. 

We provide a comparison of the three spectral models 
in Table With D and i fixed in the KERRBB and 
BHSPEC, each model has the same number of free pa- 
rameters: two for the soft/thermal component, one for 
the non-thermal component, and one for the interven- 
ing absorption column for a total of four free parame- 
ters. DISKBB provides a considerably poorer fit than 
either KERRBB or BHSPEC. The relativistically broad- 
ened spectra are a much better representation of the soft 
thermal emission than the narrower DISKBB. 

The quality of the DISKBB fit is sensitive to the 
model of the non-thermal emission. If we treat r as 
a free parameter, r drops and the fit improves slightly 
ixt = 320/176) but remains poor compared with KER- 
RBB and BHSPEC. The COMPTT spectrum steepens 
as r decreases. This extra flux in the 'tail' of the ther- 
mal component compensates for a decrease in Tin which 
allows DISKBB to better approximate the low energy 
photons. The fit further improves to xt — 2^\l\l^ if 
we replace COMPTT with a power law. This provides 
a slightly better fit than KERRBB but BHSPEC is stiU 
preferred. The best fit model now requires a steep power 
law component P 2.8 which is consistent with the 
best fit model of Haardt et al. (2001). However, the 
power law fiux now exceeds the DISKBB flux at low en- 
ergies. This result is unphysical in a picture where the 
soft X-ray emission provides the bulk of the seed pho- 
tons for the non-thermal component. This is also likely 



the explanation for why the Haardt et al. (2001) fits 
require a larger neutral hydrogen column. This incon- 
sistency was also pointed out by Yao et al. (2005) who 
find a self-consistent fit with a Comptonized multitcm- 
perature blackbody model. In contrast, the inclusion of 
non-thermal emission has little effect on the values 
for KERRBB and BHSPEC. Thus, the BeppoSAX data 
can be completely accounted for by a bare accretion disk 
spectra as long as relativistic effects on the spectra are 
included. 

BHSPEC provides a better fit (Ax^ = -32) to the data 
than KERRBB for an inclination i = 67°, a source dis- 
tance D = 52 kpc, and / = 1.7. Allowing / to vary from 
1.5-1.9 does not improve the KERRBB quality of fit sig- 
nificantly. The prescription for relativistic effects in the 
two models are essentially identical, so the differences in 
the spectral shapes are primarily due to the different pre- 
scriptions for the disk surface emission. The annuli spec- 
tra which make up the BHSPEC model have imprints 
from metal opacities and may differ from color corrected 
blackbodies by several percent. Additionally, annuli at 
different radii have local spectra which are best approx- 
imated by different values of /, with / usually being 
higher for the hotter, inner annuli. KERRBB assumes 
one value of / for the whole disk. Though these dis- 
crepancies are not at a level which is significantly greater 
than the intrinsic uncertainties in the BHSPEC model, it 
is suggestive that a model which includes atomic physics 
provides a better fit. 

Despite these differences in the quality of fit, KERRBB 
and BHSPEC both give values of a, - 0.3 for i = 67°. 
The best fit value for a, for BHSPEC is a function of a 
with a = 0.1 giving a lower a* than a = 0.01. This is the 
case in all fits and is most simply understood by exam- 
ining how changes in the parameters either harden (in- 
crease the mean photon energy) or soften (lower the mean 
photon energy) the spectra. As a* increases, the inner 
radius of the disk decreases. This results in a larger frac- 
tion of the gravitational binding energy being released 
in a smaller area of the disk surface. The resulting in- 
crease in Toff in these annuli produces a spectrum with 
higher average photon energies. Therefore, increasing a* 
hardens the spectra, even at fixed luminosity. 

The sensitivity of the spectrum to a is more complex. 
It is strongest at high M in radiation pressure domi- 
nated annuli where the surface density E is low. A larger 
a yields a lower E, making the disk less effectively opti- 
cally thick. For a range of M, the a = 0.01 annuli remain 
very effectively optically thick while the a = 0.1 annuli 
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Fig. 3.— The 66%, 90%, and 99% confidence contours in the 
a*-D plane of the best fit BHSPEC model {a = 0.01, i = 67°, 
N free) to the BeppoSAX LMC X-3 data. Here we have let the 
normalization vary by 20% above and below its nominal value of 
N = (10kpc/52kpc)2 = 0.0370. The model normaUzation N is 
enumerated on the upper horizontal axis. The vertical dashed lines 
mark uncertainty hmits associated the distance estimate in Table 

m 



become less effectively optically thick and eventually ef- 
fectively optically thin as M increases. For a = 0.1, 
the densities are lower and the temperatures are higher. 
The photons cannot thermalize as well, causing the spec- 
trum to harden significantly. At lower M, both models 
have sufhciently large E so that spectral formation occurs 
nearer the disk surface at approximately the same den- 
sities and temperatures. The a — 0.1 models still tend 
to be slightly less dense in the spectral forming region 
and are therefore slightly harder, but the differences are 
significantly smaller than at higher M. Therefore, the 
dependence of the best-fit a, on a results from an in- 
crease in a from 0.01 to 0.1 hardening the spectrum so 
that a, must be reduced to compensate. 

The best fit a* is also a function of i. As can be seen in 
Table|21 making i and D free parameters does not signifi- 
cantly improve the quality of fit in these cases. However, 
it does greatly increase the uncertainty in the best fit a, . 
This is illustrated in Figure |21 where we plot the 66%, 
90%, and 99% joint confidence contours for a, and i. 
The strong correlation exists because lowering i softens 
the spectrum so that a* must increase to compensate. 
This is partly because the line-of-sight projection of the 
azimuthal fluid velocity component decreases. This re- 
duces the blueshift and beaming of the emission from the 
approaching side of the disk, and moves the 'position' of 
the high energy tail to lower energies. Also, the projected 
disk area increases which produces a larger flux for the 
observer at infinity. This needs to be compensated by 
a decrease in M, and therefore £. This decrease in £ 
lowers Toff and again softens the spectrum. Thus, our 
ability to constrain a* is generally improved by precise, 
reliable estimates for i. The contours in Figure El suggest 
a* ~ 0.55 ± 0.2 at 90% confidence. If the binary incli- 
nation uncertainties are accurate and the inner accretion 
disk is aligned with the binary orbit, these constraints 
imply 0.2 < a, < 0.4. 

We also examine the variation of a* with the source 



distance D. The BHSPEC model normalization N is de- 
fined so that D — 10 /^/N kpc. The confidence range 
for the distance to LMC X-3 listed in Table ^ provides 
a tight constraint on N. However, N also depends on 
the absolute flux calibration of the detector so that any 
uncertainty in absolute flux translates into an effective 
uncertainty for D. Therefore, we consider fits where 
N is free to vary by 20% from its nominal values of 
N = (10kpc/52kpc)2 ^ 0.0370. The 66%, 90%, and 
99% confidence contours in the a^-D plane are shown 
in Figure O The best fit model lies at the upper limit 
of the allowed range of N. The best fit a, ~ 0.38 is 
slightly larger than the range (a* ^ 0.25 ± 0.05 ) consis- 
tent with the distance constraints which are plotted as 
vertical dashed lines in the figure. 

The strong anticorrelation between a* and D seen in 
Figure 131 exists because an increase in D leads to a de- 
crease in the flux expected at the detector (i.e. a lower 
N) . This must be accounted for by an increase in the lu- 
minosity {£) which would shift the spectral peak to higher 
energies at fixed a,. However, a, is a free parameter and 
it can be lowered so that the spectral peak remains fixed 
while £ increases. 

We also use these models to test for the possibility of 
magnetic torques on the inner accretion disk. The energy 
release by a torque increases the fraction of emission at 
small radii and increases the effective temperature of the 
annuli. This produces a hardening of the spectrum sim- 
ilar to an increase in a*. At the time of publication, 
BHSPEC spectra have only been computed from disks 
with non-zero torques for a* = 0. KERRBB can be used 
to examine torques at all a*, but the fits provide little 
constraint as Ary/77 varies over the entire range of the 
model from zero to one at 90% confidence when i is a 
free parameter. As expected, increases in Ary/ry are off- 
set by decreases in a*. An upper- limit on torque can be 
obtained by fitting BHSPEC at a* = 0. The best fit 
AT]/r] = 3 ± 0.8 with xVi^ = 245/175 and i ^ 73° ± 1°. 
This value of is only slightly greater than in the un- 
torqued case, and the best fit inclination is consistent 
with the constraints on the binary inclination, so it is 
difficult to rule out the possibility of large torques from 
these data. 

2.2.2. RXTE Data 

GD04 have already selected a sample of disk domi- 
nated RXTE observations for several sources, including 
LMC X-3. From these, we have selected a subset of 10 
epochs which evenly cover the range of disk luminosi- 
ties inferred from the GD04 analysis. The luminosities 
of these epochs are plotted versus the maximum color 
temperature in the first panel of Figure ^ The black 
filled circles correspond to the data sets used in our work 
and red triangles represent the other epochs in the GD04 
sample. These plots were generated by taking DISKBB 
fit results and making corrections for the temperature 
profile and relativistic effects (Zhang et al. 1997). This 
plot only includes epochs in which the disk component 
is inferred to account for greater than 85% of the bolo- 
metric flux. A detailed explanation of the analysis can 
be found in GD04. The values of Tdisk/^Edd and Tmax 
are evaluated using the estimates in Table so there is 
some uncertainty in collective position of these symbols 
on the plot. However, the placement of the points rela- 
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TABLE 2 
LMC X-3 BeppoSAX Fit Summary 



Model^ 


a 


i 

(deg) 


D 
(kpc) 


a* 


(keV) 


Nh 
(102"cm-2) 


xl 


DISKBB 

KERRBB 

BHSPEC 

BHSPEC 

BHSPEC 


0.1 
0.01 
0.01 


67 
67 
67 


52 
52 
52 

^-•-■^-o.o 


u.cSDciy_Q Q013 

U.l^i_0.020 
2';8+° °l^ 

U.ZOO_g Qj^g 

^4+011 
*J-3^-0.10 


±.u±oy_Q Qogo 


5.08t»i^ 

O.DO_Q 26 

0.0i5_Q 26 
02+0.31 


336/177 
275/177 
243/177 
246/177 
238/175 



Note. — All uncertainties are 90% confidence for one parameter. Parameters reported without uncertainties were held fixed during the 
fit. 

''The full Xspec model is WABS*(Model+COMPTT). 



tive to each other is robust to these uncertainties so that 
reproduction of the shapes of these L — T relations pro- 
vides an important test for our disk models. We repeat 
the same procedure for observations of J1550 and J1655 
and plot these in the center and right panels of Figure 0] 
respectively. 

As stated in section ^ the luminosity is roughly pro- 
portional to the fourth power of the maximum temper- 
ature. The dashed curves represent lines of constant / 
where L cx T^ax- I'^ the bottom panels of Figure 01 we 
have also plotted L/T^^^ in order to more easily eval- 
uate spectral hardening relative to this overall trend. 
Comparison of the data with these curves shows some 
evidence for hardening with increasing L for J1665 and 
J 1550. LMC X-3 is roughly consistent with a constant 
/, but a close examination suggests there might be weak 
signs of hardening above ~ 0.9 keV and softening at the 
highest temperatures. 

We investigate this spectral evolution with M by fit- 
ting the models directly to the data. We consider the 
same models as in section 2 .11 but we now fix the ab- 
sorption column since it is not well constrained with- 
out the low energy coverage. For LMC X-3 we fix it at 
Nh = 5.5 X 10^" cm~^ to be consistent with the Bep- 
poSAX Hts. We initially fix i and D, and fit only a single 
value of a* for all epochs. We also fix / = 1.7 for KER- 
RBB and a = 0.1 or 0.01 for BHSPEC. Only M (or 
£) is allowed to vary for each epoch. We also consider 
models with DISKBB, fitting a single normalization si- 
multaneously to all data sets. This is also consistent with 
assuming a fixed color correction and constant effective 
area for each epoch. With these choices, each model has 
the same number of free parameters. There is a single 
parameter shared by all data sets (a* or DISKBB nor- 
malization) , and two parameters for each individual data 
set: one for the soft/thermal component (M, £, or Tin) 
and a normalization for the non-thermal component. For 
BHSPEC, we also consider fits where i and D (or N) are 
free parameters. Since only a single value of either pa- 
rameter is fit for all epochs, this provides at most two 
additional parameters. 

In oder to visualize the variation of the spectral shape 
with M, we plot in Figure^ the L — T relations (solid 
curves) derived from the best fit BHSPEC models at 
fixed i for a=0.1 and 0.01. These curves are calculated 



by generating artificial spectra with our best fit mod- 
els, and then fitting them using the same procedure that 
GD04 used for the real data. The plot therefore pro- 
vides a comparison of fits with DISKBB, both to the 
data, and to artificial spectra generated from the best-fit 
BHSPEC models. It is not a direct comparison of BH- 
SPEC with the data. This explains why the best-fit BH- 
SPEC curves do not go through the 'data points' in the 
J1550 and J1655 plots. The two types of fits find a dif- 
ferent partition of the spectra between the soft/thermal 
and hard/non-thermal components, with additional fiux 
accounted for by the non-thermal emission in the BH- 
SPEC fits. The implications of this are discussed further 
in section r?.3l For comparison, we also show curves with 
-^disk/^Edd fx T^^x (dottcd lincs). These curves repre- 
sent simple, analytic estimates for the spectral harden- 
ing in effectively optically thick disks when the effects 
of Comptonization are negligible. The derivation of this 
relation in presented in the Appendix. 

Based on the discussion above, it is reasonable to ex- 
pect that the difference in quality of fit between mod- 
els with fixed / (KERRBB and DISKBB) and BHSPEC 
might be dominated by the differing predictions for the 
variation of spectral hardening with changing £. The lack 
of hardening at the highest i seen in Figure 01 for LMC 
X-3 suggests that models with constant / or BHSPEC 
with a low value of a would provide the best represen- 
tations of the data. These predictions are borne out in 
simultaneous fits to the LMC X-3 data, which are sum- 
marized in Table 13 The unfolded spectra are plotted in 
Figure El with the best fit BHSPEC model for a = 0.01 
and i = 67°. DISKBB, which is representative of a disk 
with a fixed emitting area and constant /, provides the 
best fit. A comparable is provided by KERRBB with 
i = 67° and / = 1.7. The fit with BHSPEC for a = 0.01 
gives an acceptable x^, but provides a poorer represen- 
tation than the fixed / models. The a = 0.1 model does 
not provide an acceptable fit. For BHSPEC, the largest 
contributions to come from the highest luminosity 
epoch. As seen in Figure 0] the BHSPEC models seem 
to harden too rapidly to accommodate all epochs simul- 
taneously. BHSPEC hardens more rapidly with increas- 
ing £ for a = 0.1 than 0.01, leading to the significantly 
poorer fit. In order to make the inner disk annuli more 
effectively thick, we increased the metal abundances to 
three times the solar value. However, this had limited 
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0.5 1.0 0.3 1.0 0.5 1.0 

kT,,„, (keV) kT„,,, (kcV) kT„,,, (keV) 

Fig. 4. — The disk luminosity (top panels) as a function of maximum temperature for J1655, J1550, and LMC X-3. Each symbol represents 
a DISKBB fit to one RXTE data set. These data were presented in GD04 and the reader is referred there for a complete discussion of their 
spectral analysis. The black filled circles represent the epochs which were used in this work and the red triangles represent the remaining 
GD04 data sets. The dashed curves are lines of constant color correction corresponding to / = 1.6,1.8,2.0,2.2,2.4,2.6, and 2.8 for a 
Schwarzschild black hole (see eq. 3 of GD04). The dotted lines represent curves with Ldisk/i'Edd <^ ^max- This L — T relation follows from 
simple, analytic estimates for the spectral hardening in effectively optically thick disks when the effects of Comptonization are negligible. 
The normalization is chosen arbitrarily to compare with the BHSPEC model curves. The derivation of this relation and its relevance to 
the models is discussed in the appendix. The green and blue curves show the evolution expected from the best fit BHSPEC models for i 
fixed at the binary estimate (see Table0 for a = 0.1 and 0.01, respectively. These curves were created by producing synthetic data sets 
with the BHSPEC model and replicating the spectral analysis performed by GD04. To more easily evaluate the spectral hardening relative 
to a fixed /, we have plotted L/T^^^^ (lower panels). The units on the vertical coordinates are arbitrary, but the lines of constant / (now 
horizontal) are retained for reference when comparing with the top panel. 



impact on the spectral shape and did not improve the 
quality of fit appreciably. 

The best fit values of a* and their 90% confidence inter- 
vals are summarized in Table 13 They are systematically 
lower (a* < 0.1) than the values inferred from the Bep- 
poSAX fits (a* ~ 0.3). Most of this discrepancy can be 
accounted for by cross calibration differences between the 
two observatories. Cross-calibration'^ campaigns on 3C 
273 show that the RXTE PCA flux is about 20% higher 
than the BeppoSAX data. At fixed normalization, this 
requires a 20% increase in £ which leads to an ~ 5% in- 
crease in Tcff. For a given i, this change must be offset 
by a decrease in a* which leads to a systematically lower 
value for the RXTE fits relative to BeppoSAX fits. 

2.3. XTE J1550-564 

As seen in the middle panels of Figure 01 J1550 varies 
over an order of magnitude in £ and we have chosen 
10 observations which sample this range. We fit J1550 
(and J1655) with the same model which we applied to 
LMC X-3, but we found statistically significant residuals 
consistent with reflection features. (These were appar- 
ently unnecessary in LMC X-3 because the spectra are so 
strongly disk dominated.) To account for these residuals 
we add a GAUSSIAN and apply a SMEDGE (Ebisawa 

^ http: / /heasarc. gsfc.nasa.gov/docs/asca/calibration/3c273_results 



et al. 1991) to the COMPTT component to approxi- 
mate reprocessing at the disk surface. We fix the widths 
of these components to 0.5 keV and 7 keV, respectively. 
This adds two free parameters from each new component. 
We also let the optical depth vary in the COMPTT com- 
ponent for a total of five additional parameters in each 
epoch. The fit results with these additional components 
are summarized in Table 13 

As with LMC X-3, all of the models provide an accept- 
able fit, though there is still considerable variation in the 
values. There is a statistically significant preference 
for the DISKBB model over KERRBB. Since both mod- 
els have a constant /, the difference in the quality of fit 
must be due to the overall spectral shape and not sim- 
ply its evolution with M. A comparison of the best fit 
spectral shapes shows that relativistic KERRBB mod- 
els are harder than those of DISKBB over the RXTE 
band, which seems to be primarily the result of relativis- 
tic broadening. 

The BHSPEC model with a = 0.01 provides a better 
fit than a = 0.1, consistent with the prediction of the 
L — T comparison in Figure 2| However, the a = 0.01 
model still provides a poorer fit than DISKBB. In con- 
trast to KERRBB, the best fit spectra fall off more 
steeply with increasing photon energy than the DISKBB 
spectra. This behavior seems to be primarily due to ab- 
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sorption features (primarily Fe K) in the tail of spectrum. 
When we let i and D float, the BHSPEC fit improves and 
the values are now slightly better than DISKBB. The 
best fit inclination is lower {i = 42^13 deg) and the spin is 
higher (a, — 0.72to;oi), producing shghtly harder spec- 
tra with less pronounced absorption features than in the 
i = 72° case. For i = 72°, the best fit a^, is relatively low 
(< 0.1) for both KERRBB and BHSPEC, so allowing i 
to vary changes a, considerably. 

Apparent motion > 2c has been claimed to be observed 
in the radio emission from this source (Hannikainen et 
al. 2001) suggesting i < 50° for ballistic motion. This 
implies a misalignment of at least 20° if the jet is aligned 
with the black hole angular momentum vector. The best 
fit inclination i = 42^^3 deg is consistent with this upper 
limit and may be compatible with an inner disk aligned 
with the black hole via the Bardeen-Petterson effect. 

2.4. GRO J 1655-40 

Of the three sources considered in this work, J1655 
displays the strongest evidence of hardening in Figure 
01 suggesting that the models with fixed / will provide 
poorer fits than BHSPEC. In fact, KERRBB does not 
provide an adequate fit to the data, and a luminosity 
dependent trend can be seen in the residuals. However, 
DISKBB can still provide an adequate fit and is even 
preferred to BHSPEC for i = 70°. As in J1550, the best 
fit KERRBB spectra are harder than those of DISKBB. 
These results again suggest that the overall differences 
in spectral shape (as opposed to the variation with M) 
provide the dominant effect on the quality of fit. Unlike 
LMC X-3 and J1550, the BHSPEC model with a = 0.1 
provides a better fit than a — 0.01. Comparison with 
the bottom-right panel of Figure^suggests that the soft, 
thermal component in J1655 is hardening more strongly 
with increasing luminosity than can be easily accounted 
for with the a = 0.01 model. 

At fixed inclination, BHSPEC with a = 0.1 provides 
the only relativistic fit which is marginally acceptable. 
The best fit a* = 0.62 ± 0.01 in this case suggests a 
moderate spin which is reasonably consistent with other 
investigations (a, ~ 0.7 — 0.9, Gierlinski et al. 2001; 
a, ~ 0.65 - 0.75, Shafee et al. 2006). KERRBB also 
yields a similar spin (a* = 0.6015to;oo23) fo'' / — 
The small discrepancy with the Shafee et al. (2006) is 
likely due to our choice of a single / and to differences 
in our modeling of the non-thermal emission. 

Allowing j to be a free parameter significantly reduces 
for both values of a, but a = 0.01 now provides a 
slightly better fit. As with LMC X-3 and J1550, the spin 
is very sensitive to the inclination. For both values of a 
the best fit a* = 0, the lower limit of the model. (We 
have not yet extended BHSPEC to retrograde spins.) We 
find xl = 248/302 and i = 83.8° ± 0.6° for a = 0.1, and 
Xl = 234/302 and i = 85.6° ± 0.4° for a = 0.01. 

Radio observations of J 1655 have inferred a jet inclina- 
tion of 85° ± 2° to the line of sight (Hjellming & Rupen 
1995). Assuming the jets are aligned with the angular 
momentum axis of the black hole, this would imply an 
inner disk inclination (to the plane of the sky) oii = 85°, 
consistent with the best fit i above. However, disk align- 
ment only occurs for black holes with non-zero angular 
momentum and the transition radius is expected to in- 



crease with increasing spin (Bardeen & Petterson, 1975). 
Therefore, because these fits find a, = 0, they do not 
provide a self-consistent picture for a misaligned disk sce- 
nario. 

3. DISCUSSION 

One of the principle aims of this work was to test the 
applicability of the relativistic a-disk model in BHBs. 
The spectral fitting discussed in section[21presents mixed 
results. The significant improvements in x^ relative to 
DISKBB resuhing from the BHSPEC and KERRBB fits 
to the BeppoSAX data provide a strong case for rel- 
ativistic broadening. DISKBB alone is too narrow to 
adequately approximate the soft, thermal emission from 
LMC X-3. The additional quality of fit improvement for 
BHSPEC relative to KERRBB might also be evidence 
for modified blackbody and smeared absorption features 
in the spectrum. 

In light of these results, it is surprising that DISKBB 
with a fixed normalization seems to provide a better fit to 
the RXTE data than KERRBB or BHSPEC in all three 
sources when we fix i at the binary inclination. Since 
it is preferred to both relativistic models, the difference 
cannot simply be due to the differences in the degree of 
spectral hardening as luminosity changes. This could be 
taken as evidence against relativistic broadening, but the 
innermost radii implied by the DISKBB fits are consis- 
tent with coming from near the black hole. A comparison 
of the best fit spectral shapes for all three models shows 
differences at the < 10% level. The KERRBB spectral 
shapes tend to be broader (harder in the 3-20 keV band) 
than both DISKBB and BHSPEC. The differences be- 
tween BHSPEC and KERRBB seem to be mostly due 
to broad absorption features which cause the BHSPEC 
model to fall off more strongly with increasing energy 
in the tail of the spectrum. Therefore, it is conceivable 
that DISKBB spectrum could be mimicking similar, but 
slightly weaker, features in data, though it is surprising 
that it does this consistently and effectively in all three 
sources. 

Alternatively, it may simply be that our estimates for 
D, M, or i are in error. When we allow i to be a free 
parameter, BHSPEC provides a better fit than DISKBB 
in both J1550 and J1655. In both cases, observations of 
the radio jets suggest the angular momentum of the black 
hole is misaligned with that of the binary. It is suggestive 
that in both cases we find values for i consistent with the 
constraints implied by the jets, rather than the binary 
inclination. In the case of LMC X-3, where a jet has not 
been observed, the best fit i is more nearly face on than 
the measured binary inclination. It is consistent with the 
binary inclination at 90% confidence for the BeppoSAX 
data, but not for the RXTE spectral fits. 

A third possibility is that the non-thermal emission 
and Compton reflection components are not being cor- 
rectly accounted for by our prescription. If this is a prob- 
lem, it should be minimized by looking at LMC X-3, 
which has the most disk dominated spectra of the three 
sources and relatively little evidence for reflected emis- 
sion. For LMC X-3 the BHSPEC residuals are clearly 
dominated by the most luminous epoch for which BH- 
SPEC predicts too much hardening with increasing lu- 
minosity for either value of a. If we ignore the most 
luminous epoch xt improves to 257/388 for BHSPEC 
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TABLE 3 
RXTE Fit Summary 



ModeP 


a 


a* 


LMC X-3 

xl 


XTE J1550-564 

a* xl 


GRO J1655-40 

a.* xl 


DISKBB*^ 

KERRBB<= 

BHSPEC 

BHSPEC 

BHSPEC^ 


0.1 
0.01 
0.01 


„ , -,Q+0.013 

U.113_g gj^g 


ri+0.006 

n 790+0. 036 
0. '^8_0.018 


296/431 
296/431 
1190/431 
359/431 
307/429 


097+0 005 
n+0.0055 
^-0 

„ ,,p.+0.030 

72+0 15 
'^-0.01 


230/355 
301/355 
324/355 
256/355 
221/353 


6015+0 0013 
^■°l^'-0.006 

639+° °^2 
u.uoa_g QQg 

0.00^0021 


284/304 
439/304 
330/304 
369/304 
234/302 



Note. — All uncertainties are 90% confidence for one parameter. Parameters reported without uncertainties were held fixed during the 
fit. Unless otherwise noted, the values of i, D, and M were fixed at the estimates given in Table [Tl 

''The full Xspec model is WABS*(Model+COMPTT) for LMC X-3, and WABS*{Model+GAUSSIAN+SMEDGE*COMPTT) for J1550 
and J1655. 

I'A single normalization is fit for all data sets. 

'^The hardening factor is fixed at / = 1.7. Fit parameters were selected so that limb darkening was included but self-irradiation was not. 

''Both D and i are free parameters with the value of i allowed to vary over the full range but D constrained to lie within the confidence 
limits reported in Table[Tl The best fit values are i = 43.4l^-^ deg, D = 51.4^°' kpc for LMC X-3; i = ^2^1^ deg D = 6.3510 ± 0.0080 
kpc for J1550; and i = 85.60lo:38 deg, D = 3.4000^^ gggg kpc for J1655. 



optically thin at small radii. Once the hottest annuli 
are effectively optically thin, they become increasingly 
isothermal or inverted in their temperature profiles as 
inverse Compton scattering in the now hotter surface 
layers increasingly dominates the cooling. This is much 
less efficient than thermal cooling and the spectra harden 
rapidly as temperatures rise with increasing Tcff (Davis 
& Hubeny 2006). Such effects are responsible for the 
hardening in the a = 0.1 models at high i in Figure 2| 
Our results suggest that these multi-epoch fits are sensi- 
tive to these effects and may even be able to differentiate 
between the a-disk prescription and more general models 
of angular momentum transport in these disks. 

In the context of the a prescription, our fit results seem 
to rule out fixed values of a > 0.1 for LMC X-3 and pos- 
sibly J1550. At first sight, this appears inconsistent with 
the a > 0.1 values inferred for the outburst phases of 
dwarf novae (Lasota 2001) and soft X-ray transients (e.g. 
Dubus et al. 2001). However, it is important to note 
that the disk instabilities that drive the outburst time 
scales are associated with regions of the disk where gas 
pressure dominates radiation pressure. In contrast, the 
X-ray spectra are dominated by the innermost regions of 
the disk where radiation pressure can be important at 
high M. There is no reason to believe that either a or 
the stress prescription should be the same at all radii in 
the disk. 

It may be that the classical stress prescription of equa- 
tion |(5J is not valid in radiation pressure dominated 
disks. Alternative stress prescriptions have long been 
considered, partly because they can produce thermally 
and viscously stable disks (e.g. Piran 1978). In addition, 
there have been proposals that magnetohydrodynamical 
turbulence may produce stresses that are limited to val- 
ues that are related in some way to the gas pressure (e.g. 
Sakimoto & Coroniti 1981, Merloni 2003). 

It is noteworthy that LMC X-3 reaches the highest 
Eddington ratios among the sources fit here, and that 
the spectrum even appears to soften slightly at these 
highest luminosities. A substantial reduction in stress 



with a — 0.01 and i — 67°, but only improves slightly to 
256/388 for DISKBB, providing comparable fits. Allow- 
ing i to be a free parameter allows the BHSPEC fit to 
improve even further for slightly more face on values. If i 
remains fixed at 67° , the fit with BHSPEC also improves 
by allowing the model normalization to vary within 20% 
of the nominal value. This is a larger range of normaliza- 
tion than that associated with distance uncertainty and 
accounts for possible errors in the absolute flux calibra- 
tion. Thus for LMC X-3, the shape of the spectra seem 
to be best represented by BHSPEC, but the evolution of 
the spectral hardening with £ at the highest luminosities 
is not consistent with the predictions of a simple a-disk 
model. 

3.1. Spectral Hardening and the Stress Prescription 

The three L — T relations presented in Figure 0] show 
that differences exist in the spectral evolution with disk 
luminosity from source to source. LMC X-3 has the most 
scatter and is reasonably consistent with a constant /, 
though the lower-left panel of FigureQlshows evidence of 
a weak hardening for 0.1 < £ < 0.3 and softening for £ > 
0.3. J1550 is consistent with weak hardening and J1655 
seems to show more significant hardening. A comparison 
of the KERRBB and BHSPEC fit results seems to agree 
with these descriptions. At fixed i, KERRBB with / = 
1.7 is preferred for LMC X-3, but BHSPEC provides a 
better fit in both J1550 and J1655. Fitting a single / for 
all epochs with KERRBB does not alter this result. A 
comparison of BHSPEC fits with a = 0.1 and 0.01 is also 
generally consistent. With i fixed at the estimate for the 
binary inclination, a = 0.01 is preferred for LMC X-3 
and J1550, but a = 0.1 provides a better fit for J1655. 
If I is a free parameter, a — 0.01 provides a better fit for 
all three sources, though the improvement is small for 
J1655. 

As discussed in section 12.2.11 and the Appendix, the 
sensitivity of the spectra to a comes about primarily 
because a determines E. For sufficiently large S, the 
spectral shape depends only weakly on a. However, if S 
drops sufficiently, the disks begin to become effectively 
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Energy (keV) 

Fig. 5. — The unfolded spectra for the RXTE observations of 
LMC X-3 using the best fit BHSPEC models for i = 67° , D = 52 
kpc,and a = 0.(D1. The total model component (green, solid curve), 
BHSPEC (red, long-dashed curve) and COMPTT (violet, short- 
dashed curve) are plotted. 



at ~ 0.6 — 0.7 might explain this, perhaps linked to the 
onset of a disk instability, which appears to be present 
in the only other source to consistently exceed such lu- 
minosities, GRS 1915-^105. Fits with BHSPEC to GRS 
1915-1-105 show that this source has 'stable' disk spectra 
(i.e. constant for longer than 16 sec) from I — 0.5 — 0.6 
and from t = 1 — 2 which are consistent with the expected 
hardening with luminosity (Middleton et al. 2006). How- 
ever, there are no stable disk spectra from this source in 
the range £ ~ 0.7 - 0.9, exactly where the LMC X-3 
spectra show slightly different properties than expected. 

Other effects might be important at the high lumi- 
nosities. The disk models used in BHSPEC are actually 
somewhat inconsistent for i > 0.3, as the innermost an- 
nuli then have H/R > 0.1. Radial transport of accretion 
power is increasingly important in this regime. In addi- 
tion, magnetic torques across the ISCO may be more im- 
portant (Afshordi & Paczyhski 2003). However, both of 
these effects would tend to make the inner annuli hotter 
and/or more effectively optically thin. We would expect 
this to increase the hardening, in contrast to what is ob- 
served. As the Eddington ratio increases, increased inho- 
mogeneities in the magnetorotational turbulence (Turner 
et al. 2002) and photon bubbles (Begelman 2001) may 
mitigate this by producing a softer spectrum and a ge- 
ometrically thinner disk than would be expected in a 
homogeneous model. 

3.2. Estimates for Black Hole Spins 

The total mass accreted over the lifetime of these 
sources expected to be small (King & Kolb 1999; Shafee 
et al. 2006). As a result, only a small increase in the 
angular momentum of the black hole is expected, and 
so our spin measurements are likely probing the natal 
spin distribution of the binaries. Therefore, these low to 
moderate spin estimates (a, < 0.8) place constraints on 
black hole formation scenarios. They may also constrain 



spin-dependent models for jet production (Middleton et 
al. 2006) since both J1550 and J1655 are microquasars. 
We characterize these spins as 'moderate' because even 
for a* ~ 0.8, the proximity of the ISCO to the event 
horizon and the resulting radiative efficiency {R ^ 3i?g, 
r/ ~ 0.12) are substantially less extreme than in the max- 
imally spinning spacetime (a* ^ 0.998, R ^ 1.24i?g, 
77 ~ 0.32). These moderate spins are in contrast to es- 
timates for spin at or near a* = 0.998 that have been 
inferred from other methods, including spectral fits to 
the Fe Ka line and some interpretations of the high fre- 
quency quasi-periodic oscillations (QPOs). 

Broad Fe Ka line emission has been seen in ASCA ob- 
servations of J 1550 and J 1655 (Miller et al. 2004). The 
emission was modeled with a relativistic disk line pro- 
file calculated in a maximally spinning Kerr spacetime 
(LAOR in Xspec, Laor 1991), but with inner radius al- 
lowed to vary. The best fit inner radius for J1655 was 
small with i?in < 2i?g, suggesting a* > 0.9 and possibly 
near maximal (a* ~ 0.998). J1550 is less well constrained 
with i?in ^ 4 — 6Rg depending on the model, though it 
may also be consistent with near maximal spin. The 
strongest constraints on a* come from the presence of a 
broad, asymmetric red wing which extends down to ~4 
keV in the best-fit models. Therefore, a principle source 
of uncertainty for applying this method is modeling the 
underlying continuum to accurately gauge the shape and 
extent of the line wing. We refer the reader to Done & 
Gierlinski, (2005) for a more detailed discussion of un- 
certainties associated with this method. 

The reproducibility of the frequencies of the pair of 
high frequency QPOs from one observation to the next 
suggests they might also provide a direct probe of the 
black hole spacetime. As a result, prospective models 
of QPOs often provide constraints on a*. For example, 
if the lower frequency member of the pair in J1550 and 
J 1655 is identified as an axisymmetric radial epicyclic os- 
cillation, then black hole spins of a* > 0.9 are required 
(Rezzolla et al. 2003, Torok 2005). The reason is simple: 
the radial epicyclic frequency has a maximum at some 
radius, and that maximum is below the observed QPO 
frequency for the observed black hole masses unless the 
spin is high. The same conclusion holds in diskoseismol- 
ogy models if the lower frequency member is identified 
with a low order axisymmetric "g-mode" , because it has a 
frequency less than the radial epicyclic frequency within 
the mode trapping region (e.g. Wagoner et al. 2001). 
However, there are many possible oscillation modes in 
accretion disks, and other mode identifications can be 
made which are more consistent with moderate spins 
(e.g. Blaes et al 2006a). 

In principle, the high signal-to-noise RXTE data allow 
us to place very tight constraints on the spins of BHBs 
(see e.g. Table 3). However, the relatively small uncer- 
tainties in Table 3 do not account for uncertainties in 
M, i, and D. There are degeneracies among the fitting 
parameters in the way that they affect the spectrum, e.g. 
the correlation between i and a* seen in FigureEl There- 
fore, when i and D are free parameters, a* can change 
significantly. For the RXTE data, the uncertainty ranges 
for increase when i and D are free, but still remain 
relatively small. This also leads to rather small formal 
uncertainties for i and D, though D is at the limit of 
allowed range for LMC X-3 and J1655. So these data 
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are capable of constraining all the parameters simulta- 
neously, because changes of only a few percent in the 
spectral shape in the high energy tail of the spectrum 
lead to substantial changes in x^. However, we caution 
that the models themselves are uncertain at the few per- 
cent level due to our interpolation method alone (Davis 
& Hubeny 2006). Therefore, this method of spin estima- 
tion can only be used with good precision when reliable 
and precise constraints on M, i, and D are available. 

All of these methods of spin estimation have their 
weaknesses, as all rely to varying degrees on uncertain 
physical assumptions. An important example relevant 
for both Fe Ka line and our continuum spectral fits is 
the assumption that the disk emission extends to the 
ISCO, and effectively ceases interior to this radius. In 
principle, the disk (or its emission) could be truncated 
at larger radius and our spins would be underestimates. 
Alternatively, significant emission might be generated or 
reprocessed inside the ISCO (Krolik & Hawley 2002), 
making the interpretation of both estimation methods 
more difficult. 

Several uncertainties also remain in other physical as- 
sumptions which underly the BHSPEC model, because 
we still lack a complete understanding of the magnetohy- 
drodynamical structure of the accretion flows. There are 
several modifications which likely lead to a hardening 
of the spectra, including increased dissipation near the 
disk surface (Davis et al. 2005), an increase in the den- 
sity scale height due to magnetic pressure support (Blaes 
et al. 2006b), surface irradiation by the non-thermal 
emission, and torques on the inner edge of the disk. 
Other processes might lead to a softening of the spectra. 
Opacity due to bound-bound transitions of metal ions, 
which is not included in BHSPEC should increase the 
ratio of absorption to electron scattering opacity, push- 
ing the spectrum closer to blackbody. Inhomogeneities, 
such as those caused by compressible magnetohydrody- 
namical turbulence (Turner et al. 2002) or the photon 
bubble instability (Turner et al. 2005), may also soften 
the spectrum. They may increase the effective ratio of 
absorption to scattering opacity because photon-matter 
interactions are dominated by the densest regions (Davis 
et al. 2004), or through a reduction in the density scale 
height (see e.g. Begelman 2001) leading to an increase in 
the average density of the disk interior. Given all these 
possibilities, it is difhcult to say with certainty what net 
effect modifying our assumptions or including these ad- 
ditional processes would have on the disk spectra. If 
the effects which harden the spectra are more important, 
BHSPEC will underestimate the spectral hardening and 
would then require higher spins to fit the data, making 
our spin measurements overestimates and vice-versa. 

Given these uncertainties, it is conceivable that one 
could reconcile the spectral constraints with a* ~ 0.9 
if the BHSPEC model overestimates the actual spectral 
hardening. Spins this high would bring our estimates 
in line with the published uncertainties for the Fe Ka 
estimates and some of the QPO-based measurements. 
Reconciling our spectral models with the extreme spin 
(a* ~ 0.998) is more difficult. We have attempted to 
quantify this in the case of J1655 by fitting the data 
with KERRBB for a* = 0.998, but with / as a free pa- 
rameter. For i — 70° or 85°, an adequate fit can be 
obtained only for / ~ 1 which corresponds to blackbody 



radiation. Electron scattering opacity dominates at the 
relevant temperatures so nearly blackbody emission is 
highly unlikely. One could obtain agreement with more 
reasonable color corrections (/ ~ 1.5) by allowing i to 
be a free parameter. However, this requires i < 40°, in 
disagreement with the inclinations from both the binary 
and jet observations. 

3.3. Uncertainties due to Hard X-ray Emission 

A potential difficulty for deriving constraints on a or 
a* by this method is the need to account for the non- 
thermal emission. This is particularly true with RXTE 
data when only the high energy tail of the spectrum is 
typically observed. Even though we infer the models' 
bolometric fluxes to be dominated by the softer thermal 
component, the non-thermal component accounts for a 
substantial fraction of the 3-20 keV emission in many 
of the cases. The decomposition of thermal and non- 
thermal emission is clearly dependent on the choice of 
model for the thermal emission. It can be seen in Figure 
Hthat the Tmax values derived from the best-fit BHSPEC 
models for J1550 and J1655 are softer (lower T^ax) at 
fixed idisk/^Edd than those derived by GD04 who fit 
DISKBB directly to the data. This means that emission 
which was being accounted for by the DISKBB model is 
partly being accounted for by the non-thermal emission 
in the BHSPEC fits. The fits also depends on the choice 
of model for the non-thermal component. The COMPTT 
model assumes a single Wien spectrum for the soft pho- 
ton input. A model with a multitcmperature disk spec- 
trum for the seed photon input (THCOMP, Zdziarski et 
al. 1996) provides more low energy photons for a given 
temperature. Therefore, a THCOMP spectrum which 
matches COMPTT at higher photon energies will tend 
to have more flux at lower energies. We find that the 
quality of fit can change significantly (e.g. DISKBB is 
no longer preferred to BHSPEC for J 1550 at fixed i) if 
we replace COMPTT with THCOMP, though the best- 
fit spins and preferred a values seem more robust. The 
LMC X-3 observations are much more disk dominated 
than those of J1550 and J1655, and are much less sensi- 
tive to the choice of model for the non-thermal emission. 

The effects of the non-thermal emission can be mini- 
mized if we observe these sources with detectors sensitive 
to lower photon energies, as most of the spectra peak be- 
low the lower hmit of the RXTE band. Thus, the RXTE 
fits can be sensitive to changes of only a few percent in 
both the non-thermal emission model and the shape of 
the high energy tail of the BHSPEC models. The shape 
of the high energy tail of the BHSPEC spectra is uncer- 
tain at the several percent level due to our interpolation 
scheme (Davis & Hubeny 2006) . Line-of-sight absorption 
at softer-photon energies will put practical limits on such 
modeling, and makes sources located out of the Galactic 
plane (such as LMC X-3) particularly well suited for this 
type of study. 

4. CONCLUSIONS 

We analyze disk-dominated spectra of three black hole 
binaries: LMC X-3, XTE J1550-564, and GRO J1655- 
40, fitting them with a simple multitcmperature black- 
body model (DISKBB), as well as sophisticated relativis- 
tic disk models (KERRBB and BHSPEC). For LMC X- 
3, this includes BeppoSAX AaXa., which cover the 0.1-10 
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keV energy band in which the majority of the bolometric 
flux of the disk is emitted. In this case we find a statis- 
tically significant preference for the relativistic models 
over DISKBB. At lower significance, we also find a pref- 
erence for BHSPEC over KERRBB which may suggest 
the spectra are sensitive to atomic and radiative transfer 
physics which are calculated explicitly in BHSPEC. 

We also examine RXTE spectra for each of the 
three BHBs using simultaneous, multi-epoch fits to each 
source. When we fix the relativistic model at the inde- 
pendent estimates for the binary inclination in Tabled 
we find DISKBB is preferred over both relativistic mod- 
els for both J1550 and J1655, though is still accept- 
able in most of the relativistic model fits. If we allow 
the inclination to be a free parameter, BHSPEC is the 
best-fit model for both sources. The best fit inclinations 
are both consistent with constraints inferred from radio 
observations of the jets in these sources which might be 
accounted for by a misalignment of the black hole spin 
with the binary orbital angular momentum. BHSPEC is 
also the best fit model for LMC X-3, if we ignore the high- 
est luminosity epoch where additional physics appears to 
be important. The best-fit inclination is marginally con- 
sistent with the constraints on the binary inclination in 
this source. 

Using the binary inclination estimates in Table ^ we 
are able to derive precise estimates for the black hole 
spin. However, the inferred values of the spin are func- 
tions of inclination, and the spin changes if the inclina- 
tion of the X-ray emitting, inner disk annuli differs signif- 
icantly from independent estimates. The best fit spin is 
also sensitive to the source distance and absolute flux cal- 
ibration through the model normalization. Accurate and 
precise estimates for all these parameters as well as the 
black hole mass are therefore a prerequisite for accurate 
and precise spin estimates. We find relatively moderate 
spins (a* < 0.8), even when inclination is a free param- 
eter in our fits. For J1655, our maximum spin estimate 
is only sHghtly lower than the hmits (a, > 0.9 at 90% 
confidence) implied by fits to Fe Ka lines (Miller et al. 
2004) and certain models of high frequency QPOs ( e.g. 
Wagoner et al. 2001; RezzoUa et al. 2003; Torok 2005). 
The spin of J1550 is more weakly constrained by the Fe 



Ka fits (Miller et al. 2004), and is consistent with our 
estimates. However, for both sources the best fit Fe Ka 
models are also consistent with nearly 'maximal' spins 
(a* ~ 0.998) which would be difficult to reconcile if the 
BHSPEC models provide an accurate approximation to 
the spectra of the accretion flows in these sources. 

We also find that our fits are sensitive to the assumed 
form of the angular momentum transport through its ef- 
fects on the disk surface density. We consider an a stress 
prescription with a = 0.1 and 0.01, finding approximate 
qualitative agreement between the L — T diagrams and 
the model predictions (see figure 01 . We find that all 
three BHBs are consistent with a single value of a, pre- 
ferring a — 0.01 to a = 0.1 (though only weakly in the 
case of J1655). These results are in contrast to the stan- 
dard disk instability model of soft X-ray transients (see 
e.g. Dubus et al. 2001) which requires a > 0.1 in the 
outer disk. If we include the most luminous epoch in 
our fits to the RXTE data of LMC X-3, models with 
constant / provide a better fit than BHSPEC with ci- 
ther a. BHSPEC fails to account for the most luminous 
epoch because it predicts continued spectral hardening at 
the highest Eddington ratios while the observed spectra 
appear to soften. This suggests the onset of additional 
physics which softens the spectrum as the Eddington ra- 
tio nears unity. 



We gratefully acknowledge Ivan Hubeny for his instru- 
mental contributions to the development of BHSPEC. 
We thank Eric Agol for useful discussions, and for mak- 
ing KERRTRANS publicly available. We also thank 
Marek Gierlinski, Julian Krolik, Ramesh Narayan, Re- 
becca Shafee, and Aristotle Socrates for useful discus- 
sions and/or assistance. We are grateful for comments 
on the manuscript from the anonymous referee and Jeff 
McClintock which lead to significant improvements. Part 
of this work was completed while the authors were hosted 
by the Kavli Institute for Theoretical Physics at UCSB. 
This work was supported by NASA grant NAG5-13228, 
and by the National Science Foundation under grant 
PHY99-07949. 



APPENDIX 



A SIMPLE ESTmATE FOR THE SPECTRAL HARDENING 

A major focus of this work is understanding if a standard accretion disk model can reproduce the spectra of the 
thermal dominant state, including its variation with luminosity. This task is difficult because it involves a self-consistent 
calculation of the vertical structure and radiative transfer in a number of annuli. This requires the solution of many 
coupled non-linear differential equations, and sophisticated methods (e.g. Hubeny & Lanz 1995) are needed to solve 
this problem with precision. However, it is insightful to first examine the extent to which the results of the detailed 
calculation can be inferred from simple (albeit somewhat crude) arguments before comparing with the data. 

We focus on a single annulus (the one with the largest T^s) which remains at a fixed radius R in the disk as £ varies. 
We make the further approximation that all the spectral variation is encapsulated in a single parameter /, the color 
correction in equation^] In the standard model £ — L/LEdd = / ^Edd, with ^Edd, MEdd, and rj fixed as M varies. 
We can relate M to T^s via 

where a is the Stefan-Boltzmann constant, G is the gravitational constant, and Q Keplerian frequency. For simplicity 
we have dropped the 'correction factors' due to relativistic effects and the no torque inner boundary condition. Now 
we have reduced the problem of calculating the L — T relation to finding / as a function of T^g. 
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In order to evaluate / we need an approximate scheme for the spectral formation. In our models, the spectra 
harden because electron scattering dominates the opacity, leading to deviations from local thermodynamic equilibrium 
(LTE) and modified blackbody spectra. The resulting spectra are harder than a blackbody at Toff as typical photon 
energies are higher. The photons are thermally emitted deeper in the atmosphere and their energy distributions 
are characteristic of the temperatures where they are formed. If the opacities are frequency dependent the depth of 
formation (thermalization surface) varies with frequency and photons at different frequencies originate from different 
depths with different temperatures. However, over a broad range of temperatures where bound-free processes dominate, 
the true absorption opacity is a much weaker function of frequency than in the free-free case. This 'grey' opacity 
dependence yields a spectrum which can be crudely approximated by a diluted Planck function with a temperature 
evaluated at the depth of formation r, (Davis et al. 2005). Therefore, a simple estimate for / is given by / ~ T^/Teff. 
If the opacity has a stronger frequency dependence, the Planckian shape will be a poorer approximation, but this 
estimate of / will still be crudely correct if is evaluated near the spectral peak. 

In the electron scattering dominated atmospheres we are considering, ~ TO*Kes where Kos is the approximately 
constant electron scattering opacity and m* is the column mass where the effective optical depth is equal to unity. It 
can be approximated by m, = 1/^/3k^^^k^, where Kab is the absorption opacity evaluated at m*. In general, is a 
frequency dependent quantity and should be evaluated by integrating from the surface inward. For our approximate 
estimates, it is useful to ignore the frequency dependence and consider all quantities with an asterisk subscript as 
being evaluated at a frequency near the spectral peak. For simplicity, we assume an absorption opacity of the form 
Kah = kqP^T~^ and find 

r^^^^^"\-^'/V\ (A2) 

The temperature T at an optical depth t in an annulus can be estimated using an LTE-grey model with mean 
opacities (Hubeny 1990), yielding 



■i r / ^2 



r- — +3-1/2 



Ttot / 3I]k_b 



(A3) 



Here, rtot is the optical depth at the disk midplane, and kb is the Planck mean opacity. We can use this expression 
to evaluate the temperature at the depth of formation 

^ ^T^cffT* (A4) 

where we have assumed t» ^ Ttot, t* > 3""^/^, and E 3> k^^. These limits are appropriate for annuli with large S 
and correspond to lower values of M and/or a. In this approximation, equation (|A4p shows that / depends only on 
r* through 

T. 13 



Next, we need an expression for density. In the range of Teff, E, and f2 of primary interest to us, the annuli are 
radiation pressure dominated at the midplane. However, strong gas pressure gradients are still needed to support 
the atmosphere near the surface (Hubeny 1990) where the radiation force ceases to increase with height above the 
midplane z as fast as the tidal gravity. Hydrostatic equilibrium is then given by 

cl^^n\z-zo) (A6) 

where zq is the height above which the gas pressure gradient dominates (of order the radiation pressure scale height), 
Cg = fcBT,/(/iTOp) is the local isothermal sound speed, /i is the mean molecular weight, and nip is the proton mass. 
Using the definition of column mass dm — —pdz, equation ljA6|l becomes 

^^^-{z~zo)dz. (A7) 
We can use this to solve for p, approximating the atmosphere as isothermal. We find 

p = Pocxp . (A8) 



2d 

This expression can be integrated to find m*, 

m^^Pf) 



PQC. 




(A9) 
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where Az = — and erfc denotes the compUmentary error function. The argument of t he erfc is typicaUy large in 
the range of interest. In the Umit of large x, erfc(x) ~ exp(— a::^)/(\/7rx) and equation (jA9|l can be evaluated to find 

2 

T, = Kos^Tl* — ■ (AlOj 

The expression for Az is difficult to approximate, and we must appeal to the full atmosphere calculations for guidance. 
Typically, both z^ and zq increase with Tcs but their difference tends to decrease. It decreases only weakly in annuli 
where r* < 100 and S > 10^ g cm~^, but more rapidly with Toff in annuli where t, is larger and/or S is smaller. 
However, in these annuli the assumptions underlying equation (|A10I) are increasing invalid, so over the range of interest 
we approximate Az as a constant. Including the decrease in Az would make / a stronger function of T^g, but only 

weakly since oc Az~^ and / oc r,^^^. Since E cx a~^, this approximation is better for lower a. 

Combining equations (|A2|) . HA4p . and IjAlOp . we can solve for the dependence of T*, r,, and on Teg: 
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where n = 3.5 would be the relevant value for a Kramer's opacity law. The bound-free opacity dominates the true 
absorption opacity in the models at the temperatures of interest and the mean opacity is not well approximated by a 
Kramer's law at these temperatures and densities. We find empirically that a slightly weaker temperature dependence 
{n ^ 2) yields better agreement between equation (|A2(I and the simulation results. For n = 2 we find from equation 

(jSU that / cx T^^^ oc T^^^, so that we expect / to increase with temperature, but with a weak dependence. The 
resulting L—T relation is plotted as a dotted line in Figure^ where we have used equation IjAip to obtain L/L-^dd oc T^. 
The normalization is arbitrary and chosen to facilitate comparison with the models. 

It is worth noting the role that density gradients near the surface play in producing this trend. The choice of n = 2 
yields independent of T^s, a reasonably good approximation for model annuli with large S. If we instead assume 
(as is commonly done) that radiation pressure dominated annuli are constant density with equal to the midplane 
value of the density p^md, we find a much stronger density dependence since Pmid cx (Shakura & Sunyaev 

1973). This produces a much stronger dependence of r* and / on T^s- It also yields a stronger dependence of / on 
a since Pmid a^^ . The — pmid approximation becomes increasingly relevant when S declines, either due to an 
increase in M or an increase in a (see e.g. Shakura & Sunyaev 1973). However, the large /'s that would result are 
not typically seen in the data. As pointed out by Shimura & Takahara (1995), this is largely the result of Compton 
scattering effects when and r, become large. 

It is difficult to address the effects of Compton scattering in detail without attempting a full solution of the radiative 
transfer equation. However, some insight can be gained by examining the behavior of the y-parameter (see e.g. Rybicki 
& Lightman 1979) 

y^ ^k^T-h. ^ (All) 

where T is temperature, h is Planck's constant, m,, is the mass of the electron, c is the speed of light, and N is the 
number of electron scatterings. The typical ratio of final to initial photon energy is E{/Ei ^ exp(?/). In an annulus, 
the temperature increases with optical depth as in equation (|A3|) . Photons produced near r* will tend to have higher 
energies than the electrons they encounter as they scatter out of the atmosphere. Therefore, hv > AkBT so that the 
final photon energy is less than the initial photon energy and y is negative. 

For photons produced near r,, hi/ ^ AkBT^: and N ^ for a photon reaching the photosphere. Following Hubeny 
et al. (2001), we therefore define the 'effective y-parameter' as 

y. ^ ~'^rl (A12) 

The effects of Comptonization decrease the photon energies and bring them closer to LTE near the photosphere, 
reducing the spectral hardening and lowering /. Instead of in creasing indefinitely with i ncrea sing T», / will 'saturate' 
when Ei/Ei ~ 1// or y„ ~ — In / ~ —1. Using equation (|A5|) . we can rewrite equation HA12|I as 

For ^ —1 and fceT'cff ~ 1 keV, we find / ^ 1.6, close to the commonly assume value of 1.7 (Shimura & Takahara 
1995). Above this value, / will increase only weakly because further increase in the typical energies of photons emitted 
at is offset by reductions due to the effects of Compton scattering. The precise value of / obtained will, of course, 
depend on the details, but the strong / dependence in equation (|A13|I demonstrates why the spectra of effectively 
optically thick annuli are never consistent with / ^ 2. 
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We emphasize that the arguments discussed above are only approximate with varying accuracy over the parameter 
space of interest. They are not a substitute for full atmosphere calculations. Nevertheless, comparison with these 
calculations (see e.g. Figure (S)) demonstrates that they are qualitatively correct and are useful for understanding why 
the spectral hardening depends so weakly on T^g in most BHBs. 
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